Probing Anderson localization of light via decay rate statistics 
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Abstract 

We have studied the distribution of resonant widths P (T) in one-, two- and three dimensional 
multiple light scattering systems. P (T) should follow a universal power law P (T) ~ T _1 in the 
localized regime as confirmed by extensive numerical calculations. This behavior can be interpreted 
as an unambiguous signature of exponential Anderson localization of light in open systems. 
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The research on Anderson localization of light has been of great interest 111 since it was 
originally proposed as the optical counterpart of electronic localization |2J. Localization, as 
proposed by Anderson, is defined as an inhibition of wave diffusion in infinite disordered 
media due to interference of multiple scattered waves |3[ . A much stronger definition is that 
the eigenf unctions in an infinite disordered medium are characterized by an exponential 
decay in space, \ip (r)| ~ exp (|r — r'| /£), where £ is the localization length. In finite, open 
media, waves can "leak" through the sample boundaries. Anderson localization must thus 
relate to manifestations of leakage in observables quantities. For optical systems, they are 
typically the emerging intensity, the total transmission or the coherent backscattering cone. 
The observation of an exponential scaling of transmission 0,0], as well as the rounding of the 
backscattering cone may not have definitively established localization since absorption 
could be responsible for these same effects. 

There are several criteria to determine the onset of the localized regime. The Ioffe-Regel 
criterion states that, in three dimensions (3D), localization occurs for kl ~ 1 (with k the 
light wavenumber inside the medium and i the mean free path). Another approach to define 
electromagnetic localization is based on the variance of fluctuations of transmission, even 
in the presence of action & In open s y8tem8 , the instates" ar e r_es with a 
finite energy width T (or, equivalently, with a finite lifetime t ~ due to leakage. The 
Thouless criterion establishes that localization can be set to occur when the typical time 
that an excitation needs to propagate through the entire system of size R, tx ~ 1/TV ~ 
R 2 /D (Thouless time), exceeds the maximal time scale of the system, tu ~ ~ 1/AE 

(He.enoe^tnneja.He^isthed^on — andA^e mean level spacing. 

This Thouless criterion applies to the average leakage width. Hence it is reasonable to 
assume that the statistical properties of resonance widths are strongly affected by localiza- 
tion. The aim of the present paper is to investigate how localization manifests itself in the 
distribution of resonance widths P{T) in multiple light scattering in open systems. We will 
show that P(T) exhibits the universal power law P(T) ~ T _1 in ID, 2D and 3D optical dis- 

q n 

ordered systems, thereby generalizing recent theoretical jSj and numerical |]J| studies in ID 
models of mesoscopic transport. We assert that the algebraic decay P(T) ~ represents 
a universal property of Anderson localization of light. 

Although the statistical properties of resonance widths in open systems have been exten- 



sively studied over the last years, in particular for chaotic/ballistic systems 
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12, 131, their 



2 



behavior for systems exhibiting localization has received considerably less attention. As ar- 
gued by Casati et al, P(T) should follow a power law P(T) ~ in localized, classically 
chaotic systems P(T) was analytically obtained for ID disordered systems, showing 

a slightly different power law P(T) ~ r -1,25 j^. This prediction was corroborated later by 
numerical calculations in ID and quasi-lD tight binding models (lfj. The P(T) ~ be- 
havior was also reported in ID and 3D strongly driven atomic Rydberg states in the context 
of dynamical localization ^3|. Exactly at the Anderson transition and in the diffusive 
regime P(T) was shown to follow a power law with a power different from —1. Concern- 
ing the study of P(T) for optical systems, the only work on the subject is, to the best of our 



knowledge, due to Patra 



18| , who mainly focused on the small T regime and its application 

is different from a power law, for both 
18| . It should be emphasized that the 



to random lasers. For small T, it is known that P(T 

t he d iffu8i veQQand the ,„ CallZ ed regim e S Q' 
power law decay of P(T) is expected to occur only for intermediate values of T, i.e., for not 
too large T, where the resonances are strongly coupled to the continuum and where P(T) 

[in 

decays faster than algebraically llfj, and for T > (r). 

We will present a simple physical argument, inspired by Refs. 14| . to explain the 



universal P(T) ~ behavior for the localized regime, i.e., independent on the dimen- 
sionality of the system. Due to the opening of the system, exponentially localized eigen- 
states of the corresponding closed system (linear size R) acquire a finite frequency width 
T', r" ~ e~ r '^, with r' the distance to the boundaries. Near the system edges, the leakage 
is strong and the resonances are large compared to Ft- On the other hand, far from the 
boundaries the leakage is small and the typical T in this region is much smaller than IY- 
Assuming that the resonances are uniformly distributed in space, it follows that the (inte- 
grated) probability of finding a resonance width T smaller than T', Pi nt (T < V), is equal to 
the probability of finding a resonance situated at a distance r from the boundaries larger 
than r', P (r > r'), i.e., P int (T < V) = P (r > r'). Since P (r > r') oc fi D (R - r') /fi D (R) 
with /i£> the D-dimensional volume, we conclude that the probability density is 

fiD (R - r') 



K P dr' 



fiD {R) J' 

The purely geometrical factor d/dr' (i? — r') / fin (R)} depends on the dimensionality of 
the system but does not affect the exponent in T. 



To test the validity of (JIJ in the context of Anderson localization of light, we will consider 
scalar wave propagation in disordered media using the method introduced by Rusek and 
Orlowski |l9j, |20|. This approach is based on the analysis of the Green matrices spectra, 
which describe light scattering from randomly distributed pointlike dipoles (i.e., particles 
much smaller than the wavelength of light). For an incident plane wave ipo (r) in a system 
of N identical dipoles with scattering matrix t, the field acting in the dipole at rj is given 
by 

N 

^(r i )=^ (r,.)+tE G (^)^(r i ). (2) 

The complex- valued N x N matrix G (r^) describes light propagation of the wave scattered 
by the dipole at to the dipole at rj. Since the eigenvalues Am of M = I — tG and Xg of 
G are related by Xm = 1 — t\a, and t depends on frequency u via the scattering phase shift 
8 (u) 21 1) an eigenvalue Xq with ReAc = — 1 will permit an appropriate choice for 8 (u) 
in order to obtain an eigenvalue Am = 0, which would correspond to a genuinely localized 
state somewhere inside the random medium 19]. Assuming a Breit-Wigner model for the 
scatterers, (with one resonance of width T at the position uj ) and for which 5 (u) has a 
simple form, it is possible to obtain, in a good approximation, the resonance widths F via 
Xg, r/To — 1 + ReAc j^J. In the following, we will numerically diagonalize G in ID, 2D and 
3D and calculate the distribution of resonance widths P (T) using the above approximation. 

In Fig. P (r) is calculated for ID systems composed of 100 randomly distributed 
scatterers in a linear segment for two different values of the uniform optical density p: p = 1 
and p = 10 scatterers per wavelength. In ID, all eigenstates are known to be exponentially 
localized even for weak disorder and £ is of the order of the mean free path i. P (r) exhibits 
a power law with an exponent very close to —1, in good agreement with (0). In addition, 
the exponent does not change by increasing p, i.e., by decreasing £. This demonstrates that 
the algebraic decay P (T) ~ T _1 in the localized regime is valid not only for ID models 
of mesoscopic transport P, but also for our model of wave propagation in disordered 
media. At large T, P (T) decays faster than algebraically. This can be explained by the fact 
that this region is dominated by short living resonances, typically close to the boundaries, 
for which the prediction (0) breaks down. 

Fig. (J2J shows P (r) for 2D systems of 2500 scatterers randomly distributed in a square 
for p = 1 and p = 10 scatterers per wavelength squared. According to the scaling theory of 
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localization, all eigenstates in 2D closed systems should be localized, but £ scales exponen- 



tially with 



22j. The T 1 decay of P (T) in Fig. (0) is clearly visible for both values of p 



used, with an exponent very close to —1, in excellent agreement with Remark however 
that the range of the power law is broader when p is higher. Increasing p means decreasing 
t and, consequently, decreasing £. The range of the algebraic decay P (T) ~ is expected 
to be broader as more and more states become localized. Such a behavior was also reported 
in numerical calculations within the Anderson model |l0|. For large T, P (T) decays again 
faster than algebraically as in the ID case. 

In Fig. © the 3D case is considered, where P (T) is calculated for systems composed by 
1000 point scatterers randomly distributed in a sphere (radius R) for p = 1, p = 10, p = 30 
and p = 60 scatterers per wavelength cubed. In 3D, the system is expected to undergo, 
upon varying the degree of disorder, a transition from extended states to localized states. It 
is therefore interesting to investigate if and how this transition manifests itself in P (T). As 
in the 2D case, we notice that, as p increases, the range of the algebraic decay P (T) ~ T~ a 
increases. We also remark that, as p increases, the associated exponents tend more and more 
to the value —1. The exponents, obtained by a linear fit in the range where the power law 
is present, are a ~ 0.76 for p — 1, a ~ 0.83 for p = 10, a ~ 0.95 for p = 30 and a ~ 1.1 
for p = 60. This suggests, according to (Q), the onset of the localized regime for higher 
p. In fact, the Ioffe-Regel criterion for localization [hi < 1) is estimated to be satisfied for 
p > 2tc 2 « 20 for scatterers at resonance. This condition is fulfilled for p = 30 and p = 60, 
for which a is very close to 1, showing that the system with these densities are indeed in 
the localized regime and confirming that the power law P (T) ~ can be considered a 
genuine signature of Anderson localization of light. The fact that the exponents for p = 30 
and p = 60 are not exactly equal to 1 can probably be attributed to finite-size effects. Once 
again, note that P (T) decays faster than a power law for large T. 

Let us now compare the mean width of the distribution to the inverse Thouless time 
Ft = 1/tr = QDb/R 2 , where Db is the Boltzmann diffusion constant. To estimate Tt one 
should recall that in 3D, Db = ve£*/3, where I* is the transport mean free path (which is, for 
point scatterers, equal to £) and ve the energy transport velocity, ve ~ c / [l+r dwe u / r m f] [23], 
with r dwe u = l/r the dwell time in a single scattering and r m f = £/c the mean free time. 
Applying these considerations, Tt can be written as r r /r ~ 2(£/Ry. Fig. (gj) exhibits 
P (r) for the same optical densities of Fig. (JHJ) but now with T normalized to T^. For low 
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p (p = 1 and p = 10), P (r) is peaked near Ft, showing that the system is in the diffusive 
regime. Notice that there is a non-vanishing probability to find modes that live much longer 
than tx even in the diffusive regime, the so-called "prelocalized" modes 24]. As p increases, 
we observe that P (T) is no longer centered at IV and that the probability to find a mode 
with resonance width smaller than also increases. This means that, on average, the 
modes live longer than tx- At the same time, Fig. (J3J) shows that localization manifests 
itself in P (T) not only via the broadening of the power law range but also via the fact that 
the associated exponents approach to —1. We conclude again that the P (T) ~ behavior 
is an unambiguous signature of Anderson localization of light in open media. It must be 
mentioned that the present 3D study may be relevant for recent multiple light scattering 
experiments in atomic media j^, for which modeling the scatterers by pointlike dipoles 
constitutes an excellent approximation. 

In summary, we have studied the distribution of resonance widths P (T) in ID, 2D and 3D 
multiple light scattering systems composed of randomly distributed pointlike scalar dipoles. 
We have developed a simple physical argument, based on the exponential decay of localized 
eigenfunctions, to show that P (T) should follow an universal power law P (T) ~ decay 
in the localized regime. This prediction was confirmed by extensive numerical calculations 
and demonstrates that the P (T) ~ behavior can be interpreted as an unambiguous 
signature of Anderson localization of light in open media. 
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FIG. 1: The normalized distribution of resonance widths P (T) calculated for 1000 different con- 
figurations of 100 point scatterers randomly distributed in a ID segment with two different values 
of the uniform optical density p, p = 1 (full squares) and p = 10 (open circles) scatterers per 
wavelength. The dashed line corresponds to the prediction P (T) ~ T" 1 for the localized regime 
and the solid ones are just to guide the eyes. The values of T are normalized by the resonance 
width of a single dipole Tq. 

FIG. 2: P (r) calculated for up to 50 configurations of 2500 scatterers randomly distributed in 
a square for p = 1 (full squares) and p = 10 (open circles) scatterers per wavelength squared. 
The normalization of T, as well as the significance of the solid and dashed lines, is the same as in 
Fig. <d|. 

FIG. 4: P (T) as in Fig. (J2J), but now T is normalized by the Thouless frequency, Tj-- 



FIG. 3: P (r) calculated for 100 configurations of 1000 scatterers randomly distributed in a sphere 
for p = 1 (full squares), p = 10 (open circles), p = 30 (full triangles) and p = 60 (open diamonds) 
scatterers per wavelength cubed. The normalization of T, as well as the significance of the solid 
and dashed lines, is the same of Fig. ljT)l. 
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Figure 1 Pinheiro et al. 




Figure 3 Pinheiro et al. 




